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Recently a new particle physics model called Bound Dark Matter (BDM) has been pro- 
posed [l[ in which dark matter (DM) particles are massless above a threshold energy (Ec) 
and acquire mass below it due to nonperturbative methods. Therefore, the BDM model 
describes DM particles which are relativistic, hot dark matter (HDM) in the inner regions of 
galaxies and describes nonrclativistic, cold dark matter (CDM) where halo density is below 



T^ • Pc = E^. To realize this idea in galaxies we use a particular DM cored profile that contains 



three parameters: a typical scale length (r^) and density (po) of the halo, and a core radius 



O " (re) stemming from the relativistic nature of the BDM model. We test this model by fit- 

-4— > ■ ting rotation curves of seventeen Low Surface Brightness (LSB) galaxies from The HI Nearby 

j^ I Galaxy Survey (THINGS). Since the energy Ec parameterizes the phase transition due to the 

underlying particle physics model, it is independent on the details of galaxy and/or structure 

fvq . formation and therefore the DM profile parameters rs,rc,Ec are constrained, leaving only 

^ I two free parameters. The high spatial and velocity resolution of this sample allows to derive 

^^ • the model parameters through the numerical implementation of the x^-goodness-of-fit test 

ly-v . to the mass models. We compare the fittings with those of Navarro- Frenk- White (NFW), 

l/^ ■ Burkert, and Pseudo-Isothermal (ISO) profiles. Through the results we conclude that the 

r~^ ■ BDM profile fits better, or equally well, than NFW, Burkert, and ISO profiles and agree with 

^^ , previous results implying that cored profiles arc preferred over the N-body motivated cuspy 

profiles. We also compute 2D likelihoods of the BDM parameters Tc and Ec for the different 

galaxies and matter contents, and find an average galaxy core radius Tc = 300 pc and a 

transition energy between hot and cold dark matter at E,, = O.ll^QQy eV when the DM halo 

^1^ ■ is the only component, therefore the maximum dark matter contribution in galaxies. In a 

M . more realistic analysis, as in Kroupa mass model, we obtain a core re = 1.48 kpc, and energy 

■ i?. = 0.06t°:°^eV. 
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I. INTRODUCTION 

In the last decades intense work to understand the distribution of dark matter (DM) in galaxies 
has been published [jtSI]. Among the considered, late-type Low Surface Brightness (LSB) galaxies 
are of special interest since it is believed they are dominated by DM, and high resolution HI, 
Ha and optical data can help to distinguish among the different DM profiles proposed in the 
literature. There are essentially two types of profiles, the ones stemming from cosmological A^-body 
simulations that have a cusp in its inner region, e.g. Navarro- Frenk- White (NFW) profile [a, lZI]. 
On the other hand, the phenomenological motived cored profiles, such as the Burkert or Pseudo- 
Isothermal (ISO) profiles [3]. Cuspy and cored profiles can both be fitted to most LSB rotation 
curves, but with a marked preference for a cored inner region with constant density. Furthermore, 
cuspy profiles that do fit to galaxies in many cases suffer from a parameter inconsistency, since their 
concentrations (c) are too low and velocities V200 are too high in comparison to the ones expected 
from cosmological simulations [5|, l9| . Some other fits indicate that cuspy profiles require to tune 



fine the observer's line of sight with axisymmetric potentials 10|]. However, different systematics 



may play an important role in the observations such as noncircular motions, resolution of data 



and other issues llMla]. And there are attempts to reconcile both approaches through evolution 



of DM halo profiles including baryonic processes 17|, |l8| to transform cuspy to shallower profiles 



that follow the solid-body velocity curve (y ~ r) observed in late-type LSB. However, fthis issue is 



a matter of recent debate, see e.g. 19|, |20|] for a recent discussion. 



Alternatively to the above models, recently one of us proposed a profile based on a particle 
physics model called bound dark matter (BDM) [l|], in which DM particles at high densities, as in 
galactic inner regions, are relativistic, i.e. Hot DM (HDM), but in the outer regions they behave as 
standard CDM. To realize this idea in galaxies we use in the present work a particular DM cored 
profile that contains three parameters: a typical scale length (r^) and density (po) of the halo 
and a core radius [vc) stemming from the relativistic nature of the BDM model. The galactic core 
density is given by pc = E^ and the profile properties will determine the energy scale of the particle 
physics model. We will show that the LSB rotation curves yield a phase transition energy scale 
Ec, between HDM and CDM for our BDM profile, at Ec = 0.111q'o7 eV, when we consider only 
DM and Ec = O.OGlQ'gg when considering gas, DM halo, and the minimum contribution of stars 



21I ]. The Ec parameter is a new fundamental scale for DM which can be theoretically determined 
using gauge group dynamics, i.e. it does not depend on the properties of galaxies, once the gauge 
group is known. However, even though we propose Ec as a new fundamental constant for DM, 
close related to the mass of our DM particle, and it is important to note that its value is not yet 
known and we require observational evidence to determine it. The same is true for all masses of 
the standard model of particles, i.e. their value is not predicted by the standard model and it is 
the experimental results that fixed them in a consistent manner. We use here the information on 
galactic rotation curves to determine Ec- However, the extracting of Ec does depend not only on 
the quality of the observational data, the mass models used but also on the choice of DM profile. 

To perform this task, we use The HI Nearby Galaxy Survey (THINGS), which collects high 
spectral resolution data revealing extended measurements of gas rotation velocities and circular 



baryonic matter trajectories [2^. Given these properties it is adequate to test the above-mentioned 
DM profiles with THINGS, which has been used to test different core/cusp mass profiles. For disk- 
dominated galaxies the core and cusp profiles fit equally well, however for LSB galaxies there is a 



clear preference for core profiles over the cuspy models jlql. Aiialysis of different high resolution 



datasets have confirmed this tendency in past recent years [9|.I17|.|23h25|. In this work, we present a 
study of the rotation curves using BDM, NFW, Burkert and ISO DM profiles, taking into account 
the contribution from different mass models: i) DM alone; ii) DM and gas; and iii) DM, gas and 
the stellar disk. By fitting the models to the data we find that the BDM profile fits equally well or 
better than the cored profiles and only for a few LSB galaxies our model resembles that of NFW. 

The kinematics of stars bring a very challenging problem in the analysis, mainly due to the 
uncertainty of the mass-to-light ratio (T*) and to its dominant behavior close to the galactic center. 



Some considerations have been made in order to reduce this uncertainty in the parameters [26l429l| , 
but still the stellar contribution is not well known. In our case, we study four different T^, models 
and the most inner part of the data enable us to compute the BDM parameters re and Ec, both 
related to the galactic core. We present the la and 2a likelihood contour plots of these parameters 
for the different galaxies and mass models. We notice that being Vc different for each mass model it 
is consistent within the 2a error for each galaxy. Additionally, we find that the galaxy core radius 
is the order of re ~ 300 pc and the Ec ~ 0(0.1 eV) which gives the transition between hot and cold 



DM, this is shown in preliminary study [2ll |. It is also interesting to note the coincidence in the 
magnitude of the sum of neutrino masses with the magnitude of Ec obtained in this analysis, and 
how this could open an interesting connection between the generation of DM and neutrinos masses. 



but this will be presented elsewhere 30(]. The magnitude of the energy transition Ec is of the same 
order as the upper limit of the total mass of the neutrinos J2 ''^f < 0-58 eV(95% CL) from the 
seven year Wilkinson Microwave Anisotropic Probe (WMAP), Baryon Acoustic Oscillation (BAO), 



and the Hubble constant {Hq) J3ll.l32l|. In addition, the recent analysis of the Atacama Cosmology 



Telescope [33] that combined CMB data with measurements of BAO and the Hubble constant 



reported an excess of the effective number of relativistic degrees of freedom, N^ff = 4.6 it 0.8, in 



consistency with Ref. [34], opening the possibility to have extra nonstandard- model relativistic 
degrees of freedom, as our BDM particle. 

We organized this work as follows: In Sec. [II] we explain the particle model behind the BDM 
profile. In Sec. IIIII we present our BDM profile, which has NFW as a limit at low energies. In 
Sec. IIVI we describe the galaxy sample considered for the rotation curves study. The different mass 
models and components (gas, stars, DM halo) are presented in detail in Sec. [Vj while the different 
hypothesis of the mass-to-light ratio are discussed in subsection IV B[ The results and conclusions 
are presented in Sec. IVII and IVIIl respectively. For convenience, some figures and tables are 
shown in the Appendices. We show in Appendices [A] and [B] the fitting values corresponding to 
the free parameters for some of the cored profiles for the maximum stellar contribution. Finally, 
in Appendix [C] we present the rotation curves of the different profiles and galaxies and we include 
the confidence contours levels, pointing out the considerations made for each galaxy. 



II. PARTICLE MODEL 

We now present the physics and motivation behind our BDM model [l|]. Cosmological evolution 
of g auge groups, similar to QCD, have been studied to understand the nature of dark energy 



34J I and also DM 35|. A particle mass can be generated either by the Higgs mechanism or by a 



nonperturbative effect. In the Standard Model (SM) the fundamental particles (quarks, electrons, 



neutrinos) get their mass by the interaction with the Higgs field that acquires a nonvanishing 
vacuum expectation value at the electroweak scale E'ew = 0(100 GeV), while at higher energies all 
SM particles have vanishing masses. On the other hand, the nonperturbative gauge mechanism is 
based on the strength of gauge interaction and the mass of the particles is expected to be at the 
same order of magnitude as the phase transition scale as for example protons and neutrons. For 
asymptotically free gauge groups, such as the strong Quantum Chromodynamics (QCD) force in 
the SM, the gauge coupling constant becomes strong at low energies and binds the elementary par- 
ticles (quarks) forming neutral particles such as protons, neutrons, and mesons. The condensation 
or phase transition scale is defined as the energy where the gauge coupling constant g becomes 
strong, 5(A) ^ 1, giving a condensation scale Ac = Aj e~^'" '^i , where h is the one-loop beta 
function which depends only on the number of fields in the gauge group. For example for a SUSY 
gauge group SU{Nc),Nf, where Nc{Nf) is the number of colors (flavors), we have b = 3Nc — Nf 
and gi is the value of the coupling constant at an initial scale Aj. The fact that Ac is exponentially 
suppressed compared to Aj allows as to understand why Ac can be much smaller then the initial 
Ei which may be identified with the Planck, Inflation or Unification scale. In our case the gauge 
group and elementary fields are not part of the standard model (SM). Our dark gauge group is 
assumed to interact with the SM only through gravity and is widely predicted by extensions of the 
SM, such as brane or string theories. We can relate E^ to Ac since the energy density depends on 
the average energy per particle and the particle number density n, i.e. pc = E^ = Ac n. 

For asymptotic free gauge groups, such as QCD, the low energy states consist of bound gauge 
singlets, such as baryons and mesons in QCD, formed by fundamental (nearly massless) particles, 
quarks in QCD. The order of magnitude of the mass of these particles is [l|] 

mBS = dEc (1) 

with d = 0(1) a proportionality constant. In QCD one has Ec — 200 MeV with the pion mass m.,^ ~ 
140 MeV while the baryons mass (protons and neutrons) m^ ~ 940 MeV, i.e. the proportionality 
constant is in the range 0.7 < d < 5, and with bound mass much larger than the mass of the quarks 
(niu ~ (1 — 3) MeV, m^ ~ (3.5 — 6) MeV). Clearly the mass of the bound states is not the sum 
of its elementary particles but is due to the nonperturbative effects of the strong force and is well 
parameterized by Ec- The dynamical formation of bound states is not completely understood since 



it involves nonperturbative physics. However, it has been shown in RHIC [3a| that at high density, 
above the transition scale Ec, the QCD quarks do indeed behave as free particles, while at low 
energies there are no free elementary quarks and all quarks form gauge neutral bound states. Since 
the interaction strength increases at lower energies, the formation of bound states is expected to be 
larger at the smallest possible particle bound state energy Ebs (i-e. Ebs = itT'Bs) with momentum 



2 _ rp2 ^2 



E^g — rriBs — 0- '^^^ energy distribution of bound states formation is still under investigation 
37l | and for simplicity we take here p = which gives a vanishing particle velocity for the bound 
states. 

It is precisely this nonperturbative gauge mechanism that we have in mind for our bound states 
dark matter BDM. Of course, in our case the gauge group and elementary fields are not part of 
the SM. Our "dark" gauge group is assumed to interact with the SM only through gravity and is 
widely predicted by extensions of the SM, such as brane or string theories. Even though we have 
motivated our BDM in terms of a well motivated particle physics model we stress the fact that 



the cosmological implications of BDM do not depend on its origin. The BDM is defined by a DM 
that at high energy densities, Pbdrn > Pc the particles behave as relativistic HDM with a particle 
velocity v = c while for lower energy densities Pbdm. < Pc, the BDM are cold bound state particles, 
i.e. CDM with t; < c. 

There are two natural places where one may encounter high energy densities pbdm for dark matter. 
One is at early cosmological times and the second place is at the inner regions of galaxies. In the 
first case, we define Oc as the transition scale factor where 

Pc = E^ = Pbdm{a = ac) (2) 

and the Universe is dense and for a < ac the BDM particles are relativistic and redshift as radiation 
while for a > Uc we have nonrelativistic BDM particles and they redshift as matter, i.e. 

Pc < Pbdm{a < ac) oc a"^, pc > Pbdm{a > ac) oc a~^ (3) 

In the second case, away from the center of the galaxies the energy density is smaller and it increases 
towards its center and in a NFW profile it blows up p — ?> oo when the radius r — >■ 0. Therefore we 
will have for pbdm < Pc and the BDM particles are cold, i.e. CDM, but in the inner region once 
Pbdm ^ Pc one has a transition and the BDM particles become relativistic, 

Pc < Pbdmir < Tc) with V Pt, c, pc> Pbdm{r > Tc) with V ^0 (4) 

where re is defined by Pbdmif = ^c) = Pc- The dispersion velocity v of the BDM particles is vanishing 
for r larger than Vc, the particles are CDM, while at the galactic center the BDM particles are 
relativistic and for r < re it is f = c it gives a core inner region at r < Vc and with Pcore = Pc, the 
energy density of the galaxy at Vc- 



III. BDM PROFILE 

The BDM model simply consists of particles that at high energy densities are massless relativistic 
particles with the speed of light, but at low energy densities they acquire a large mass, due to non- 
perturbative quantum field effects, and become nonrelativistic with a vanishing (small) dispersion 
velocity. The phase transition takes place at an energy density defined as pc = Ec and is related 
to the mass of the DM particle. The value of pc can be determined theoretically, given a gauge 
group model, or phenomenologically by consistency with cosmological or astrophysical data. In 
the present work we estimate its value through the study of the BDM from rotation curves. It is 
worth pointing out that the mass of the SM particles are also free parameters in the SM and it is 
the physics in colliders, such as in the LHC at CERN, which defines its value. In our case, since 
the BDM particles do interact only weakly with the SM we cannot use collider data to measure its 
mass and we are left with astronomical and cosmological observations. 

The average energy densities in galactic halos is of the order pg '^ lO^pcr (Pcd being the critical 
Universe's background density) and as long as pg < pc we expect a standard CDM galaxy profile, 
which may be given by the NFW profile, pniw The NFW profile has a cuspy inner region with 
Pnfw diverging in the center of the galaxy. Therefore, once one approaches the center of the galaxy 
the energy density increases in the NFW profile and once it reaches the point pg = pc we encounter 



the phase transition and the BDM particles become massless. Inside r < r^. the BDM particles are 
relativistic and the DM energy density Tc remains constant avoiding a galactic cusp. Of course we 
would expect a smooth transition region between these two distinct behaviors but we expect the 
effect of the thickness of this transition region to be small and we will not consider it here. 

Since our BDM behaves as CDM for p < pc as long as the density of the galaxy is pg < pc we 
expect to have a NFW type profile in this limiting case. Therefore, the proposed BDM profile [l| 
is given by a cored CDM profile 

Pbdm— -. :—, r2 V>) 



fe+^)(l+^)" 

Ipc 



(6) 



with Tc < Tg and Vg^pQ are typical NFW halo dependent parameters. The BDM profile coincides 
with pnfw at large radius but has a core inner region, when the galaxy energy density phdm reaches 
the value pc = -E^ at r ~ Tc and for r^ <C r^ we have 

Pcore = PhdmV = re) = Pc = "^ , ('} 

with the core radius given by Tc = ^^p- The value of re depends on the galaxy profile parameters 
po and rg. 

If we assume that the transition energy of the BDM particles is a fundamental parameter in the 
DM nature that can be constrained by fitting the rotation curve of galaxies, we can substitute one of 
the three free parameters by the density pe in the BDM profile, Eq. ([5]), using the relation ([7]). We 
call "fixed-BDM" to the profile that is obtained as a result of substituting the free parameter pQ by 
the fundamental quantity pc as in Eq.Q. The fixed-BDM profile depends now on the fundamental 
density pc, which is the same for all galaxies, and on the two free parameters r^ and re that depend 
on the morphology of each galaxy. 

The slope of BDM profile is 

^ d log/Q _ r/re{l + 3r/rs + 2re/rs) , , 

dlogr (1 + r/rc)(l + r/r^) 

and takes the values a = (0, 1/2, 1) for r = (0, Vc, re <^ r <^ r^) and a = (2, 3) for r = {rs,rs <S r). 
If we are interested in the inner region of the galaxy, for values of r <^ r^, Eq. ([8]) can be 
approximated for Tc <C r^ with 

< a = ^— < 1 (9) 

1 + y ' 

in terms of y = r/re with < y < oo. Given a value of a we can determine the value of 

y = r/re = «/(! - a). 

For future reference (see Sec. IVip we introduce here a phenomenological ansatz for a model 
independent DM profile, useful to investigate the central region of a galaxy, given by pa = Po r~°' 
with po,a constant parameters. 



Fittings to observational data have shown that the NFW halo profile is not a g ood description 



for rotation curves, and it is gene rally preferred a core dominated halo model [2j|, although some 



scatter in a could be possible [33] • In the present work our objective is to test the realization of the 
BDM model, Eq. ([5]), through rotation curves of a wide sample of galaxies. We aim to derive the 
kinematics properties of the galaxies and, in particular, we constrain the slope a and the transition 
energy Ec- 



IV. THE SAMPLE 

With the current HI data provided by THINGS high resolution and excellent sensitivity of the 
velocity fields and rotation curves are available making possible to revisit some of the outstanding 
questions about DM, and to compare the applicability of both core or cusp models. THINGS 
galaxies have an observing data sample of 34 nearby galaxies containing a large range of luminosities 
and Hubble types at the desired sub-kpc resolution, but we limit our sample to seventeen low 
luminous (early type and dwarf) galaxies with smooth, symmetric and extended to large radii 
rotation curves and small or none bulge, see Table HI These properties provide a good estimate of 
the DM halo in galaxies because it is believed that it dominates over all other components at all 
radii. For technical details and systematic effects treatment of the observations of the THINGS 



sample refer to Refs. 22. l39l. |40|| and for a complete analysis of its rotation curves see la], hereafter 
deBlokOS. 

The mass models (in Sec. |V]) are constructed with the rotation curves extracted with THINGS 



data and the 3.6 /im data from SINGS (Spitzer Infrared Nearby Galaxies Survey) [4ll ]. We follow 
the analysis of deBlokOS and McGaugh et al. (2007) [4j] for the sample considered here. 



V. MASS MODELS 

Our mass models include the three main components of a spiral galaxy: thin gaseous disk, Vq, 
a thick stellar disk, 14, and a DM halo, Vh- In most cases the stellar disk can be well described 
by a single exponential disk. When necessary, in a small number of galaxies we have considered 
an additional central component, a bulge, containing a small fraction of the total luminosity of the 



galaxy, as described by deBlokOS la] ■ The gravitational potential of the galaxy is the sum of each 



mass component, thus the observed rotation velocity is 



Vobs' = Vh' + Vg' + T,V,'. (10) 

As an input we need the observational rotation curve, Vobs-, the rotation curve of the gas compo- 
nent, Vq, and the stellar component 14- With this data we can extract information for the DM halo 
through our theoretical model. We describe in more detail the treatment for the gas component in 
Sec. IV Al as well for the stellar component in Sec. IV B[ In Sec. IV CI we describe the DM models 
used. 



THINGS 









DIET-SALPETER KROUPA 


Galaxy 


Distance 


Rd 


logjQ Mstar logio Mstar 


(1) 


(2) 


(3) 


(4) (5) 



NGC 925 


9.2 


3.30 


10.01 


9.86 


NGC 2366 


3.4 


1.76 


8.41 


8.26 


NGC 2403 


3.2 


1.81 


9.71 


9.56 


NGC 2403 


3.2 


1.81 


9.67 


9.52 


NGC 2841 


14.1 


4.22 


11.04 


10.88 


NGC 2903 


8.9 


2.40 


10.15 


10 


NGC 2976 


3.6 


0.91 


9.25 


9.1 


NGC 3031 


3.6 


1.93 


10.84 


10.69 


NGC 3198 


13.8 


3.18 


10.4 


10.25 


NGC 3198 


13.8 


3.06 


10.45 


10.3 


NGC 3521 


10.7 


3.09 


11.09 


10.94 


NGC 3621 


6.6 


2.61 


10.29 


10.14 


NGC 4736 


4.7 


1.99 


10.27 


10.12 


NGC 5055 


10.1 


3.68 


11.09 


10.94 


NGC 6946 


5.9 


2.97 


10.77 


10.62 


NGC 7331 


14.7 


2.41 


11.22 


11.07 


NGC 7793 


3.9 


1.25 


9.44 


9.29 


IC 2574 


4 


2.56 


9.02 


8.87 


DDO 154 


4.3 


0.72 


7.42 


7.27 



TABLE I: Sample of late-type and dwarf galaxies taken from THINGS as presented in Walter et al [22l |. 
Columns: (1) Galaxy's name. (2) Distance to the galaxy in Mpc. (3) Characteristic radius of the stellar 
disk in kpc as given in deBlokOS. (4) Logarithm of the stellar mass disk when considering the diet-Salpeter 
IMF in solar masses {Mq), and (5) Logarithm of the stellar mass disk when considering the Kroupa IMF 

(Mq). 



A. Neutral Gas Distribution 

For the gas we assume an infinitely thin disk in order to compute the corresponding rotation 
curve. For more technical details we refer to deBlokOS. We point out that the case of a disk with 
sufficient central depression in the mass distribution can yield a net force pointing outwards, and 
this generates an imaginary rotation velocity and therefore a negative Vq. An imaginary velocity 
is just a reflect of the effective force of a test particle caused by a nonspherical mass distribution 
with a depression mass in the center. We have not include the contribution of the molecular gas 
since its surface density is only a few percent of the that of the stars, therefore its contribution is 
reflected in a small increase in T^, 



43|. 



B. Stellar Distribution 



To model the stellar disk we use the approximation of a radial exponential profile of zero thickness, 



the Freeman disk J44l |. since the disk vertical scale height does not change appreciably with radius 



and the correction to the velocity is around 5% in most cases [45[. Thus, the central surface 
density is S(i?) = Sqc"^'^'^, where Rd is the scale length of the disk and Sq is the central surface 
density with units [MqPC~'^]. These two parameters are obtained first by fitting the observed 
surface brightness profile, extracted from the SINGS images at the 3.6/um band and synthesized by 
deBlokOS, to the linear formula /u(i?) = ^o + 1.0857i?/i?d where ^o is the central surface brightness 
given in observational units [mag arcsec"^], //q and the surface brightness are related by a simple 



change of units. We get the surface density thanks to the mass-to- Hght ratio T^, an additional free 
parameter in the mass model, introduced because we generally can only measure the distribution 
of the light instead of the mass. 



The rotation velocity of an exponential disk is given by the well know Freeman formula 44 1 
K(2/') a y'^ [Io{y')Ko{y') - h{y')Ki{y')], where y' = R/{2Rd) and I„ and K^ are the modified 
Bessel functions of order n of the first and second kind, respectively. 

When we deduce the DM properties, the T^ as well as the Initial Mass Function (IMF) remain 
as the major sources of uncertainty. The unknown value of T^ makes it hard to pick out between 
cusp and core profiles. In order to constrain the wide range of DM halo parameters we follow 



the approach used by deBlokOS instead of the van Albada & Sancisi 26[] . We considered a "diet" 



Salpeter IMF [27f| in which stellar mass population syntheses have proved j46l | to maximize the disk 



mass contribution (maximum disk) for a given photometric constraint, J-K color band. We also 



consider the Kroupa IMF [28[] based on stellar population studies in the Milky Way that produces 
lower disk masses that minimizes the baryonic contribution by a factor of 1.4 less massive than the 
diet-Salpeter ones. 

In the numerical analysis we use the T^ from the 3.6/im obtained in deBlokOS using an empirical 
approach into the relation between the 3.6/xm emission and T^ using the J-K colors given in the 



2MASS Large Galaxy Atlas [47 1 



Disk galaxies sometimes show radial color gradients, and it is believed that this provides an 
indication of stellar population between the inner and outer regions of a galaxy and produce T^ 



gradients between these two regions of the disk J48l |. We take the T^ as a function of the radius in 



order to consider the different stars contribution as it depends on the region that we were analyzing. 



The T^ has been modeled (e.g. Salpeter [27|], Kroupa |28|], Bottema[29|]), but the precise value 
for an individual galaxy is not well known and depends on extinction, star formation history, IMF, 
among others. Some assumptions have to be made respect to T^ in order to reduce the number of 
free parameters in the model. We present a disk-halo decomposition using different assumptions 
for T^ for the galaxy sample considered here. 

Minimal disk: This model assumes that the observed rotation curve is only due to DM. This 
gives the upper limit on how concentrated the DM component can be in the galaxy. 

Minimal disk plus gas: The contribution of the atomic gas and the DM halo is taken into 
account, but stars have no contribution (T^, = 0). 

Kroupa: Studies of the stellar population in the Milky Way suggest that the Kroupa IMF produces 
low disk masses that we consider as the minimal limit for the stellar disk. 

diet-Salpeter: Here T^ is set to be a constant value based on the diet-Salpeter IMF, in which 
the stellar population synthesis model has proven to give a maximum stellar disk for a given 
photometric constraint. 

Free T^: Here we ignore a priori any knowledge of the IMF and treat T^ as an extra free parameter 
in the model, and we let the program to choose the best-fitted value for T*. 



10 



The contribution of the atomic gas is considered in the Kroupa, diet-Salpeter, and Free T^ mass 
models. In our fit we consider radial color gradient and, when present, the bulge for the stellar 
disk. In the results section we make emphasis on some difference in the core fit when we do not 
take into account one or both of these two ingredients in the stellar disk. 



C. DM Halos 

We use three well-known models for the DM distribution. ACDM simulations come up with a 
DM density profile independent of the mass of the halo characterized by a cusp central density, 
the NFW profile. On the other hand observational determinations of the inner mass density 
distribution seems to indicate that mass density profiles of DM halos can be better described using 
an approximately constant-density inner core [p ~ r" with a <C 1). This core has a typical size 



of order of a kpc [49|, l50|, and examples of these profiles are the ISO and Burkert. Here we also 



consider our proposed BDM to test it and compare it. 

Strictly speaking, we are dealing with circular rotation velocities of test particles in the plane of 
the galaxy and we assume spherical halos. For this distribution of matter the circular velocity at 
radius r is given by V^{r) = GM{r)/r. We do not consider the adiabatic contraction of the dark 
matter halo but be are aware that it predicts the increase of dark matter density in the 5% of the 



virial radius J5ll ]. 



1. NFW Halo 



The NFW profile takes the form 



^0 /'11^ 

Pnfw — — ; -2) (.J-J-j 






where r<j is the characteristic radius of the halo, and po is related to the density of the Universe at 
the time of collapse of the DM halo. This mass distribution gives rise to a halo rotation curve 



V^{r) = ^"^"'^° 



r + rs 



+ ln''" + "^ 



(12) 



This density profile has an inner and outer slope of -1 and -3, respectively. The inner slope implies 
a density cusp. The halo density can be specified in terms of a concentration parameter c = r2oo/^s 
that indicates the amount of collapse that the DM halo has undergone, where the radius r2oo is 
defined as the radii at which the density contrast of the galaxy is 200 times greater than the critical 
density pcr, defined as pcr = SH"^ /{SirG), i7 being the Hubble parameter. 

We use the po and r^ as the free parameters for the NFW model instead of the classical c and V200 j 
because these two parameters can be interpreted and compared directly with the free parameters 
of other DM halo profiles. 
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2. BDM Halo 

Recalling that the BDM model predicts an inner galaxy core radius determined by the energy 
scale Ec^ which corresponds to the phase transition energy scale of the subjacent elementary particle 
model. The proposed profile is shown in Eq. ([5]) coincides with pjif„ profile at large radii since 
p oc r~'^ but has a core inner region at r = re where p oc porsr~^. Its circular velocity is given by 



Vnirf = ^^^^ 



, + ^ # In 1 + - + 



In ( 1 + - 



(13) 

where the additional free parameter re is the galaxy core radius that demarcates the place where 
the BDM particles have an energy greater than Ec and behaves as HDM for r < re- If the radius 
is greater than Tc the BDM particles acquire a large mass through a nonperturbative mechanism 
and behave as CDM. r^ is just a characteristic scale for the DM halo. 



3. Burkert Halo 



The cored Burkert halo profile is given by 

PB 



Po 



(l + f5)[l + (^)'] 



(14) 



where po is the central density, Vg is the core radius. The rotation curves caused by this halo are 
given by 



y^(.)2 = ^^^ (in 



r. ) \ r? 



2 arctan 



(15) 



4. Pseudo-Isothermal Halo 



The spherical pseudo-isothermal (ISO) halo has a density profile 

Po 



Piso 



1+ f 



2' 



(16) 



where po here is the central density of the halo and r^ is the core radius of the halo. The corre- 
sponding DM rotation curve is 



Vnirf = AnGpor^, 



1 arctan 1 — 



(17) 



D. Computing the Mass Models 



We use the observed rotation curve, stellar, and gas component as an input for the numerical 
code, in order to obtain the properties of the DM halo. In order to fit the observational velocity 
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curve with the theoretical model we use the x^"goodness-of-fit test (x^-test), that tell us how 
"close" are the theoretical to the observed values. In general the x^-test statistics is of the form: 

^ - h V -i ) ' ^ ^ 

where a is the standard deviation, and n is the number of observations. 

We fit the free parameters of the DM halo for the Kroupa minimal disk and diet-Salpeter max- 
imum disk with a x^'test. When we minimize x^ we use different methods of minimization, 
Differential Evolution, NelderMead, and SimulatedAnnealing, in order to be sure that we are not 
in a local minimal, and at the same time we put constraints on the free parameter in order to have 
values greater than or equal to zero to obtain physical reasonable values. For BDM we constrain 
Tc to values smaller than r<j, because it makes no sense to consider core radius greater than the 
characteristic radius of the DM halo. 

Comparison of the fits derived can tell us which of the DM models is preferred. More important 
are the differences between the reduced Xred ~ X^/('^ ~ P ~ 1) values, where n is the number of 
observations and p is the number of fitted parameters. 

The uncertainties in the rotation velocity are reflected in the uncertainties in the model param- 
eters. 

As the ISO and Burkert halo have no basis in standard cosmology, there are no a priori expec- 
tations for its model parameters. One can derive them and check for possible trends with other 



fundamental galaxy parameters, as other authors have done [52 1 



VI. RESULTS 

We have analyzed the rotation curves of seventeen galaxies using five disk models (Sec. IV Bp and 
four different DM profiles: our BDM, NFW, Burkert, and pseudo-isothermal (Sec. IV Cp . Since we 
have a large number of tables and figures, and because we want to emphasize the main results, in 
this section we mainly present three disk models (min.disk, min.disk+gas, Kroupa) and two DM 
profiles (BDM and NFW). However, each DM profile is analyzed with all the disk models except 
for the free T^, for which only the BDM profile is tested. The rest of the tables and figures are in 
the Appendices as explained below. In the following we shall use the term "mass models" to refer 
to different disk mass scenarios for the analysis of the rotation curves as described in Sec. IV Bl and 
denote "min.disk" and "min.disk+gas" as minimal disk and minimal disk plus gas abbreviations, 
respectively. It is worth pointing out that including color gradients we obtain a better fit to the 
rotation curve in our analysis. We have also consider the effect of a bulge when it is present in 
the galaxy and found that the rotation curves are still consistent with a core r^ ^ 0. However, 
since the bulge contributes significantly to the inner part of the galaxy the value of Tc is smaller 
when a bulge is taken into account and in some galaxies the 1<t standard deviation of r^ may also 
consistent with a vanishing core. The adequate surface brightness distribution for the bulge has 



been proven to be of an exponential form [16j] which is the model that we consider here. Galaxies 



with a bulge deserve a closer analysis to the parameter estimation, as we discuss below. 
We present the results of NFW and BDM profiles in order to compare them since the BDM 
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TABLE II: We present the results discussed in Sec. I VII for BDM (Eq. ((Sj) with its three free parameters {re, 
Vs, and Po) and the NFW (Eq. Hill) ) profiles with the minimal disk model. Also, in the last three columns, 
we present the result obtained with the BDM profiles fixing the value of Ec to 0.11 eV as explained in Sec. 
I VI CI Column (2) shows the ratio between the core radius Vc and the scale distance rs, and galaxies are 
grouped by de value of this quotient. Columns (3-6) are the fitted results when only DM and the BDM 
profile are considered in the mass model. (7-9) show the results for NFW with the same mass model. In 
(10-12) we show the results for the fixed BDM model. All distance scales such as re and rs are given in kpc. 
Logarithm base 10 of the densities po is given in Mo/kpc^. (6,9,12) show the value of x^ is normalized to 
the numbers of data points minus the number of free parameter. 



profile that has NFW as a hmit when Tc ^^ 0. The fitted values of the profile parameters are shown 
in the following mass models order: min.disk, min.disk-|-gas, and Kroupa are in Tables H H IIIH and 
IIVI respectively. We present the la and 2a likelihood contours plots of po and Tc in Figures [2l [Sj 
H and El 

We leave in Appendix El the results for the diet-Salpeter (Table |X|), the BDM free T^ scenario 
(Table IXIIip . and the fitted values for the Burkert and ISO profiles with different mass models 
(Tables IXI V| IXVp . The details of the analysis for each galaxy as well as the figures when considering 
the different mass scenarios are in Appendix ICl 

We have grouped the galaxies into three blocks according to the fitted values of Tc of the BDM 
profile as shown in column (2) of Table [III The first one, Group A (G.A.), is composed of galaxies 
with a core radius r^ > r^ i^ 0. The second, Group B (G.B.), is formed by galaxies in which 
Tc took its upper limit value, Tc = r^, and finally, the Group C (G.C.) is where Tc is negligible, 
'^d^s < 10~^. For the sake of simplicity we conserve the same group structure for the forthcoming 
tables. We shall explain the physical interpretation of the results for each group in the following 
paragraphs. 
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TABLE III: Here we present the fitted values for the minimal disk plus gas mass model when we considered 
the pbdm and the pntv profile and the complete set of data as discussed in Sec. IVII columns (2-6) and (7-9) 
respectively. Columns (10-12) present the fitted values when fixing the energy E^ = 0.08 eV, cf lVI Al Units 

and the set of galaxies are as show in Tablellll 



Let US firstly discuss G.A. galaxies. This group comprises galaxies with reasonable fitted values 
for the BDM parameters rc,rs, and po, and it is composed of the galaxies DDO 154, NGC 2841, 
NGC 3031, NGC 3621, NGC 4736, NGC 6946, and NGC 7793. In the galaxies NGC 3031, NGC 
4736, and NGC 6946 a bulge is present and it gives a better fit, but the x^ is only slightly reduced. 
On the other hand, the bulge in galaxy NGC 3031 with a diet-Salpeter mass model is overestimated 
and the fit is worst when is analyzed without the bulge. 

For G.A. galaxies we obtain values for re < 200 pc and an average r^ ~ 4 kpc in the minimal 
disk analysis that are typical values for a galaxy, except for DDO 154 which has a larger core, 
Vc = 1.35 kpc. On the other hand, the average of the scale radius for the NFW profile is slightly 
bigger, r^ ~ 6 kpc, shifting the cusp to larger radii. In Tables H H IIIH and IIVI we show the fittings 
for BDM and NFW profiles with the three mass models mentioned above. We see that the BDM 
core radius decreases, when more mass components are taking into account. The reason for this 
is that Tc and the amount of disk and bulge masses are degenerated parameters, and this happens 
more when e.g. the stellar disk has a more dominant behavior close to the center of the galaxy, 
preventing us from having a better precision for re- 
in Appendix IC C.AI we show the fittings for all DM profiles and mass models and the rotation 
curves for each galaxy. From the likelihoods (e.g. Fig. ([2])) we see that the core radius decreases, 
and even Tc — > becomes valid, when more mass components are taking into account. The reason 
for this is that Vc and T^ are degenerated parameters, and this happens more when the stellar 
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TABLE IV: This table show the fitted values when considering the Kroupa IMF for the value of the stellar 
disk, columns (2-6) show the BDM and columns (7-9) NFW parameters. Columns (10-12) are the parameters 
obtained with the BDM profile by fixing Ec = 0.05 eV, cf. IVI Al . Units and the set of galaxies are as show 

in Table nil 



disk has a dominant behavior close to the center of the galaxy, preventing us from having a better 
precision. 

We notice that for all the G.A. galaxies the x^ in BDM is smaller than NFW's. The extra 
parameter {vc) in BDM, ho^wever, makes the reduced Xred equivalent for both profiles. The Xrcd ™ 
the BDM profile tends to NFW's value -when Tc — )• as explained in Sec. IIIII When the value of 
Tc > 40 pc, BDM is clearly better fitted than NFW, and also than the other two cored profiles, 
cf. Tables HIl IXIVI and IXVI for the min.disk mass model. Burkert and ISO profiles have difficulties 
when fitting a couple galaxies (NGC 4736, NGC 3621) having xl^d > ^ '^'^'^^^ ^DM and NFW each 
have Xred — 2- -^°^ ^^^ gcilaxy DDO 154 all cored profiles fit equally well for every mass model, 
since they have a xfcd — -*-' meaning that the fits lie inside the uncertainty of the observations. For 
galaxies NGC 3031 and NGC 7793 none of the DM profiles provide a good fit, being BDM the best 
fit for both galaxies with a Xred > ^ ^'^^ each of the mass models. Such a large Xred ^^ ^^^ ^^ ^^^ 
uncommon velocity drop in the rotation curve at large radii in these galaxies. 

We now turn to the analysis of confidence levels of re and po in the BDM profile stemming from 
the different mass models for this set of galaxies, shown in Fig. ([2]). From this table we deduce 
that the core average is Vc — 40 pc and a core energy E^ — 0.22 eV for the min.disk mass model, 
where we have taken a logarithm-normal distribution.lt is worth noticing that energy Ec has a 
smaller dispersion than re or r^ and we supports the hypothesis that a signal that Ec is a constant 
parameter. This indicates the fundamental importance of the energy parameter in the BDM model. 
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The likelihood contour plots also reveal how the central contribution of the gas makes r^ to increase 
(c.f. min.disk vs. min.disk + gas), but on the contrary, stellar disk fades away the evidence of the 
core making r^ — > 0, e.g. the galaxies NGC 3031, NGC 3621, and NGC 6946 for the Kroupa mass 
model are consistent with a null re at la. We also observe that even when re is different for each 
mass model its value lies inside the confidence levels obtained from the other mass models for the 
same galaxy. For example, for the galaxy NGC 3621 the central value of Tc < 60 pc for min.disk 
and min.disk+gas, but the confidence levels are still equivalent with a zero core. On the other 
hand, when stars are considered in Kroupa the central value is Tc -^ 0, but confidence contours 
show that can be well accepted values up to 30 pc. This, on the other hand, means that the amount 
of data in this galaxy is not good enough to clearly discern between a core or cusp profile. We 
obtain Tc = r^ for the galaxy NGC 4736 (Kroupa mass model), in which the likelihoods restrict 
quite strongly the Tc value. 

Let us now discuss Group B (G.B.) which is composed of the galaxies IC 2574, NGC 2366, NGC 
2903, NGC 2976, NGC 3198, NGC 3521, and NGC 925. These galaxies have preferred a fitted 
values Tc — Tg, that makes our BDM profile of a particular type, p oc {rc + r)~'^. The Tc value for IC 
2574, NGC 2976, and NGC 925 galaxies is of the order of the farthest observation from the galactic 
center which entails that the density of the DM halo is constant at all observed radii. In G.B. the 
only galaxies that have a bulge are NGC 2903 and NGC 3198. The galaxy NGC 3198 presents a 
sudden reddening in the inner most region that might indicate the presence of a central component, 
therefore we performed the fits with and without a bulge for this galaxy. The bulge component in 
both galaxies overestimated the rotation velocity, indicating that is likely to be dynamically less 
important than suggested by the mass models. The fits for the G.B. galaxies are in Tables H H IIIH 
IIV|, and|Xl and the rotation curves are in Appendix IC C.Bl 

It is worth noticing that for this set of galaxies the BDM profile fits much better than NFW's, 
having the former smaller Xred- ^^^ fact, Burkert and ISO profiles also fit better than NFW's, 
indicating that a core is needed for G.B. galaxies. The BDM profile fits better or equally well 
than the other two cored profiles in all the analyzed scenarios. The BDM and Burkert profiles are 
slightly better than ISO's profile when we consider min.disk and min.disk+gas analysis. Otherwise, 
for Kroupa and diet-Salpeter mass models the BDM and ISO profiles have equivalent x^d ^^^'^ ^^ 
a bit better than Burkert 's. 

The la and 2a likelihood contour plots of Tc and po for the G.B. are shown in Fig. ([3]). We notice 
that when more mass components are included in the analysis the contours plots become broader. 
The BDM profile consistent fits values of re — r^ for all mass models have an average Tc value of 6 
and 5 kpc for min.disk and min.disk+gas, respectively. There is an exception when considering the 
Kroupa and diet-Salpeter mass models in IC 2574, NGC 2903, and NGC 3521 galaxies for which 
Tc 7^ Ts , as seen in Table |IV] and |Xl 

In general, the rotation curves fits prefer cored profiles over NFW's, but for this group of galaxies 
we found that re — rg, that defines the particular cored profile mentioned before. To understand 
in more detail the parameter estimation in these galaxies we performed an analysis of the inner 
galactic region, see Sec. IVIAI The galaxies IC 2574, NGC 2976, and NGC 925 are characterized 
for having a rotation curve with constant positive slope at all radii, see Figs. [H [U and El It is 
important to emphasize that for these galaxies the different results obtained from the BDM and 
inner analysis are all in good agreement. First, the fit of BDM profile results in re — rg that is of 



17 

the order of the maximum observed radius. Then, for the inner profile pin we obtain rg/rf. > 10^ in 
most of the cases. Finally, we obtained a < 0.9 for these galaxies when fitting pa = Por~°^ taking 
into account all the observations, see Table IVIi The different results obtained from the different 
analysis imply that the DM halo density for these galaxies is constant up to the maximum observed 
radius. This clearly indicates that DM halos need to have cored profiles for the correct description 
of these galaxies. 

The la and 2(T confidence levels of re and pc are shown in Fig. ([9]). With the inner analysis 
we notice that the length of the core radius is reduced, for example, the average Tc = 2.2 kpc, 
Tc = 3.5 kpc for min.disk and min.disk+gas, respectively, with energy Ec ~ 0.06 eV in both cases 
and assuming a logarithm-normal distribution in our sample. We notice since stars are difficult to 
treat, for the Kroupa and diet-Salpeter mass models, the confidence levels increase up to an order 
of magnitude for the Tc value. 

The last and smallest group is Group C (G.C.) composed of galaxies NGC 2403, NGC 5055, and 
NGC 7331. In these galaxies is present a bulge in their stellar structure, however, their contribution 
is not enough to adequately explain the rotation curves. This is probably an indication that the 
bulge has been dynamically underestimated in the calculations, as suggested by the Kroupa or 
diet-Salpeter mass models, except for the NGC 2403 galaxy that has a good fit considering a 
diet-Salpeter mass model. These galaxies have fitted values such that Tc/ts < 10~^, c.f. Table 
im The reduced number of data close to galactic center prevents us from finding whether these 
galaxies actually have or not a core. ISO and Burkert are poorly fitting profiles, meanwhile BDM 
has a fitted value for re -^ 0, so, BDM is equivalent to NFW profile and both have equivalent 
fits. Their values are displayed in Table HH for the minimal disk mass model and rotation curves, 
and confidence contours of Tc and po are shown in Appendix [C] in Fig. @. We conclude that 
NFW/BDM profiles fit better than Burkert's and ISO's profile. 

All the mass models are consistent with a null core, except for the diet-Salpeter analysis where 
the core take values of Tc > r^ which has no physical interpretation in our model and indicates an 
overestimation of the stellar disk. Galaxies NGC 2403 and NGC 5055 present only few observational 
points close to the center, therefore we decided to perform an inner analysis that is described in 
section IVI A[ We compute the confidence level contours of Tc and pc, see Fig. (j4]). We obtain 
< Tc < 50 pc as a plausible interval that can fit the observations within the 2cr error and 
Tc < lO^^r^. Although a cusp profile is preferred for these galaxies, a more definitive conclusion 
would demand more data close to the center in these galaxies. 

To summarize, we have seen that in G.A. and G.B. cored profiles fit better than the cuspy 
NFW profile. The BDM profile prefers values for re ^ and fits better or equally well than the 
other two cored profiles. In the cases where re is very small {re < 10~^rs) a larger number of 
observational data is needed close to the center in order to clearly discern between cored or cuspy 
profiles. BDM, Burkert, and ISO profiles have good fits for G.B. galaxies, where NFW over-predicts 
the velocity in the inner parts of the galaxy or does fit with unphysical values for r^ > O(IO^) 
or pq < C'(lO^)M0/kpc'^ . In G.C. galaxies the Burkert and ISO profiles do not have a good 
fit, and the BDM profile demands a very small core radius {re < 10~^rs), but the number of 
inner data is poor and it does not give enough information on re- Assuming these data, BDM 
results are indistinguishable from NFW. Thus, BDM-NFW have much better agreement than ISO 
and Burkert profiles. Confidence level contours show consistency in the fitted values of re for the 
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different mass models. Given the richer structure the BDM profile has, through its transition from 
CDM to HDM, it allows for a better explanation of the rotation curves of the different galaxies. 

The free mass model treats T^ as an extra free parameter in the model, and the best value for 
T* is chosen by fitting the theory to the observations. The mass of the stellar disk is directly 
related to T^^ with the Freeman equation, refer to Sec. IV B[ The results when analyzing the BDM 
profile in the free mass model are presented in Table IXIIII We notice from Tables |I] and IXIIII that 
values of the stellar mass in the free mass model are in good agreement with the photometric disk 
masses derived from the 3.6/im images in combination with stellar population synthesis arguments 
predicting diet-Salpeter as the maximum stellar disk consistent with deBlokOS [16j, nevertheless 
the diet-Salpeter mass scenario overestimated the rotation velocity of most of the galaxies and 
therefore being less realistic than the Kroupa mass scenario. 
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TABLE V: We show the values obtained from the inner analysis, for all the mass model (except diet-Salpeter, 
cf. Appendix |X|, with the profile pi„. Rm is the maximum radius to consider when making the inner-core 
analysis. The core radius re is given in [kpc] and the logarithm of the central density pc in [Mo/kpc^], both 
parameters obtained from the inner analysis in which we fit the profile pin ~ 2pc{l + r/rc)^^ to the data 
below Rm where we obtained the goodness-of-fit x?od- Using the parameters (re, r^, and po) obtained from 
the analysis of the complete set of data and the pbdm, we compute xt which is the contribution of the data 

corresponding to the inner region, r < Rm- 



In this section we perform an analysis of the inner region of the galaxies. We only do it to galaxies 
in which we found a re — r^, i.e. in Group G.B. These galaxies have a large number of points in 
the outer region and the DM profile is therefore fixed by these points. However, the estimation of 
the DM profile parameters gives a bad fit to the inner region of these galaxies and gives then no 
information about the core. We will see that if we use only the galactic inner region we get a much 
better fit. 

To extract information of the galactic central region and determine the core radius, central density 
with a inner slope of the profiles, and we examine the central region data in two different scenarios: 
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TABLE VI: Fits with the profile p = por'" c.f. Eq. (|20|l for galaxies belonging to group G.B. for all the mass 
models (except diet-Salpeter, cf. Appendix |X]). The distance until which we consider the observations is 
given by Rm given in Table FVl The magnitude of the slope of the rotation curve is given by the dimensionless 
parameter a. The re is the core radius, given in Kpc, obtain based on the value of a and Eq. (|2ip. 



i) The first one by taking the limit r <C r^ in the BDM profile which turns out to be 
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ii) The second is with the ansatz profile 
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where a is a constant slope. Both approaches are related through Eq. ([9]). If the slope takes values 
a < 1/2 it implies that one is in the region where r < re, hence we have a core behavior. 

In case (i), Eq. ()19p . when the BDM profile is approximated by /9j„, we take the data points for 
r < Minimum{rc, r^) and determine the parameters r^ Pc, and the corresponding chi-square, Xinn- 
The values obtained using this method are shown in Table |Vl The Xt it is computed using also 
only the inner points but with the Vc and pc values of the complete data. We notice that x^nn i^ 
much smaller than Xt in most of the cases. 

We present the different G.B. galaxies that where analyzed when the DM halo is the only mass 
component with the inner approximation for the BDM profile pm up to a radius r < Minimumrc, Vg 
in Fig. ([5]-[8]) with a purple, long, dashed line. From this first approach we obtain the core distance 
Tc and the central density pc- The gray points are the observations with its error bars. Below we 
show the difference between the observational and the theoretical approach, the purple region 
represent the error bars and the line the fitted curve. Fitted values are presented in Table |V] and 
the confidence levels between pc and r^ in Fig. [9l 

In case (ii), Eq. (j20p . we analyze the rotation curves with pa and we fit the free parameters a and 
po using the same dataset taken for the inner analysis described in case (i). The values obtained 
for a are in good agreement, cf. Eq. ([9]), with those established for r^. in case i), see Table IVll 
For most of the galaxies we obtain values in the interval y < 1.08, that implies a dominant core 
behavior. For galaxies NGC 2976 and NGC 3198 we obtain y ~ 5.67 which means that their inner 
data correspond to distances beyond the core region, such that re < r < Vg, which can be still 
accepted for the core analysis since we include the core and the transition region between Vc and 
Ts- Given these results, we can now trust to the Vc values obtained from the pin analysis in case 
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i), since we confirm that the set of data points used for inner analysis correctly describes the core 
region of the galaxy. 

It is interesting to notice that we can obtain the value of re from the value of a. We developed 
a different approach in order to calculate the core radii from the value of the slope. The constant 
a obtained from the fit with p^ is the average value of the a{r), Eq. Q, in the interval from rmin 
to rmax, where rmin is the first observation and rmax is the maximum radius where we are doing 
the inner analysis. The average value of a is given by 



a 



rc + r 



■ dr = 1 + 



In 



Tc + rn 



Tc + Tr, 



(21) 



We can numerically solve the last equation in order to obtain the value of re- The values obtained 
from these method are shown in Table IVH from which we can see that the core radius agree in the 
order of magnitude with those obtain from the /jj„ analysis for most of the galaxies. 
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TABLE VII: We show the main parameters, Tc, pc and Ec, for the BDM model as defined in Sec. IIIII for 
the min.disk, min.disk+gas, and Kroupa mass model. For the G.A. galaxy we present the results obtained 
for analysis considering the BDM profile, Eq. ((5]). The values for G.B. galaxies are those obtained from the 
inner analysis using pi„ as described in the Sec lVIl The core r^ is given in kpc. The logarithm base 10 of 
the density pc is given in pc in MQ/kpc^. Finally, the energy Ec in which the DM transition take place is 

given in eV. 



B. BDM Energy Phase Transition E^ 



In the inner region of the galaxies and for r <^ rg the relevant parameters are po and Tc and we 
present their la and 2cr likelihood contours for the different mass models in Figures [21 [3l [H and 
m Using the fit of the rotation curves for the seventeen LSB galaxies for our BDM halo profile, 
Eq. ([5|) presented in Sec lVIl we determine the properties of the DM halo through its three free 
parameters {re, r^, and po) and constrain the energy value Ec, in which the BDM particles acquire 
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FIG. 1: We plot the values of re vs iJc for each galaxy and different mass models. Circles represent galaxies 
with re <^ rs, diamond markers are galaxies corresponding re = ra group, but considering re and Ee 
obtained from the inner analysis in Sec I VI A] In all cases the areas in yellow and blue represent the la and 

2a confidence levels. 



a non-perturbative mass and has the transition from CDM to HDM, as explained in Sec. IIIII We 
calculate the maximum, minimum, and central value of Ec and r^ for the four mass models and 
present the results in Table rvTIIl The values of the core radius re, central energy density pc and 
transition energy Ec = pj , as defined in Sec lIIH are shown in Tables rvTII and IXIII and we plot Ec 
vs Tc in Figs {Tj for the different mass models and galaxies of G.A. and G.B., taking into account 
the inner analysis of Sec iVI A|l and its respectively la and 2a likelihood contours errors. The 
circles and diamonds represent the galaxies from G.A. and G.B., respectively. We use a log-normal 
distribution for Ec,rc since these quantities are constraint from (0, oo), and the la, 2a represent a 
A X 10 standard deviation from the central value with A = Ec or A = Vc- 

The average core radius re and the energy Ec of the transition between HDM and CDM depends 
on the choice of mass models used. We can see that in general the core radius is larger for the 
min.disk+gas than for the min.disk but decreases when bulge and stars are considered as in Kroupa 
and diet-Salpeter. From Figs (|10l - ll6p we see that the contribution of the gas is very small at small 
or negligible radius, however Kroupa, and even more diet-Salpeter, have a considerable contribution 
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BDM Statistics 
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0.04 0.11 0.33 
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3.02 


0.01 0.30 7.02 


1.37 


Min.Disk+gas 


0.10 


0.05 0.08 0.15 


0.25 


8.35 


3.81 


0.15 0.98 6.60 


0.83 


Kroupa 


0.08 


0.03 0.06 0.13 


0.33 


7.65 


6.59 


0.16 1.48 13.73 


0.97 


diet-Salpeter 


0.13 


0.02 0.07 0.21 


0.47 


7.92 


5.85 


0.17 1.70 17.37 


1.01 



TABLE VIII: We show the statistics of the transition energy Ec and the core re for the BDM profile for the 
different mass models. In column (2) we show the arithmetic mean energy transition from HDM to CDM, 
Sc„. In columns (3-5) we present the median energy {Ec), and the its related uncertainty {Ec± = SclO ), 
where the a is the standard deviation of the \ogEc distribution. In columns (8-12) we present the same 
statistical parameters but related with the core distance Tc- The energy Ec and the core distance are given 
in unit of eV and kpc, respectively. In column (7) we have pc = Ec in units of Mo/kpc'^. 



of mass at small radius 0(0.1 - 1 kpc) and therefore the required contribution of DM in this range 
of radius would be smaller than in min.disk or min.disk-l-gas mass models. Therefore, it is not 
surprising that the extraction of the properties of any DM profile for small radius in Kroupa or 
diet-Salpeter gives a less constraint to r^ i.e. a larger standard deviation. 

For galaxies and mass models where the central value of Tc is vanishing, the energy density pc or 
Ec can not be calculated because it blows up as seen in Eq.d?]). In these cases we have a blank (i.e. 
" — ") in Tables IVTIl and [XIII However, even though the central value of re < Vg 10^^ it is consistent 
with a non-vanishing value at a la and 2a deviation. More data and better understanding about 
the distribution of stars in these galaxies are needed to settled the question of a core DM profile. 

The average values for Tc and E^ and the standard deviations a are presented for the different 
mass models in Table IVnll For the min.disk mass scenario we have 



£;c = 0.11 X l0±0-^6 gy ^ 0.1 ll[]:^^ eV; re 



300 X 10=^^-^'^ pc. 



(22) 



i.e. 13 pc < Tc < 5308 pc, while in the more realistic approach, the energy and the core for the 
Kroupa mass model take the values 



Ec = 0.06 X l0±°-33 eV = 0.061^-^1 eV; r^ = 1.48 x l0±°-97 kpc, 



(23) 



The figures of Ec vs re for the different mass models are shown in Fig [TJ We note that when more 
components are considered, the la and 2a confidence levels of Vc become larger and the transition 
energy Ec decreases but are all consistent within a 1 and 2a levels as seen in Tabl JVTIIl 

Clearly the scale of Ec = 0(0.1 eV) is very close the sum of neutrino masses. Therefore, an 
interesting connection could be further explored between the phase transition scale Ec of our BDM 
model, which also sets the mass of the BDM particles, to the neutrino mass generation mechanism 
30|. However, we would like to emphasize that our BDM are not neutrinos, since neutrinos are 
HDM, and a combination of CDM plus neutrinos would have a cuspy NFW profile with CDM 
dominating in the inner region of the galaxies, assuming that CDM do indeed have a NFW profile 
as suggested by N-body simulations ja, lZ|. A main difference between neutrinos and BDM is that 
neutrinos were in thermal equilibrium at S ^ MeV with the SM particles, but BDM may not have 
been in thermal equilibrium with SM at all. 
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C. Fixed-BDM profile 



Now we are able to extract further information of our BDM model. The BDM model predicts 
that the transition scale Ec is an universal constant which depends on the nature of the interaction 
of the DM particles J35l | and as explained in Sec JIIII and is independent of the galaxies details and 
therefore it should be the same for all galaxies. However, the determination on the value of Ec 
depends on the quality of the data, the choice of matter content and also on the choice of the 
galactic profile as the one proposed in Eq.([5]). 

Our BDM profile has three parameters rs,rc and pc = E^. In Sec lVIBl we presented the analysis 
and conclusions of comparing the BDM profile to the rotation curves as described in Sec. IVII We 
now proceed as follow: we obtain an average value of Ec for each of the different mass models and 
then we take this average value and proceed to reanalyze the rotation curves for the BDM profile 
but with only two free parameters, i.e. we take pc fixed in Eq.(l6|) and we have re and r^ as free 
parameters. Following, we present these results and conclusions. 

The results are shown in Tables HH- IIVI and in Table [IV] for each mass model, and Table HXl show 
the statics values of r^ and re for this analysis. We generally find that in most of the galaxies the 
fixed-BDM profile is better than NFW regardless of the mass model used. The fixed-BDM profile 
has a smaller xL ~ 1 than NFW, e.g. in DDO 154, NGC 2366, NGC 3198, or because the fit is 
clearly competitive with that of NFW, Ax^ed < !> e.g. NGC 4736, NGC 2841, NGC 6946. It is 
also relevant that the core obtained in the min.disk scenario with the fixed-BDM profile (460 pc) 
is of the order to the one obtained with the BDM profile (300 pc) 







Fixed-BDM Statistics 
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Min.Disk 


0.11 


0.81 0.11 0.46 


1.85 


0.61 


4.71 


0.51 2.64 13.44 


0.71 


Min.Disk+gas 


0.08 


1.86 0.40 1.28 


4.11 


0.51 


3.25 


1.03 2.41 5.62 


0.37 


Kroupa 


0.06 


2.01 0.44 1.27 


3.65 


0.46 


13.84 


2.90 8.56 25.26 


0.47 


diet-Salpetcr 


0.07 


1.59 0.25 0.85 


3.50 


0.62 


11.99 


2.80 7.84 22.01 


0.45 



TABLE IX: We show the statistics for the core re and r^ resulting from the fit of the fixed-BDM profile, 
Eq. ((Gjl for the difi'erent mass models and the computed Ec shown in Table IVIIII We had considered the 
log-normal distribution. In columns (3) we show the arithmetic mean fc. In columns (4-6) we present the 
median energy (re), and the its related uncertainty {rc± = rclO "), where the a is the standard deviation 
of the log re distribution. In columns (8-12) we present the same statistical parameters but related with the 

distance r^. The distance are given in unit of kpc. 



VII. CONCLUSIONS 



We have presented the analysis of the rotation curves using four different DM profiles (BDM, 
NFW, Burkert, and ISO) and five different disk mass components (only DM (min.disk), DM plus 
gas (min.disk-|-gas), DM plus gas plus stellar disk (Kroupa, diet-Salpeter, Free T^) for a subsample 
of the THINGS galaxies as described in Sec. |Vl 
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The new tested profile, BDM, is introduced in Sec. IIIII and corresponds to CDM particles below 
certain critical energy density pc = E^ while it behaves as HDM above pc with a core radius Vc- 
Therefore, the resulting profile behaves as NFW's for large radius, but it has a core in the central 
region. 

We have found that fits to the THINGS galaxies favor core profiles over the cuspy NFW in 
accordance to general results found in other works [5|, ll6| for the same galaxies. We studied the 
cases when the five different disk mass components are taking into account and conclude that BDM 
is in good agreement with observations for all the disk mass scenarios and fits better or equally well 
than Burkert and ISO profiles. For three galaxies (Group C) the fits to the BDM profile demanded 
''c/'^s < 10~^ and in these cases the x^ ^ire almost identical to NFW's. For these galaxies a cuspy 
profile fits best because there is not enough observational data close to the galactic center to test 
it. 

The observational resolution of the THINGS sample is of high quality, but one still needs data 
closer to the center of some galaxies, on scales smaller than 200 pc, in order to discern between 
cored or cuspy profiles. This is because stars pose a very challenging problem when testing the 
core-cusp problem, largely due to the uncertainty of the mass-to-light ratio and their dominant 
behavior close to the center of the galaxy. We have performed an analysis of the inner regions 
of Group B and C galaxies, for which Tc < r^, extracting information from the rotation curves 
consistent with a core profile. Confidence level contours show that the core radius depends on the 
number of disk mass components taking into account in the analysis. We found that re is highly 
constrained if the stellar disk has a dominant behavior close to the center of the galaxy. 

We computed the 1 and 2a confidence levels for the BDM parameters Tc and Ec for the different 
galaxies and mass models. We found that the energy of transition between HDM and CDM takes 
places for most galaxies with an energy Ec = 0(0.1 eV) and most of the galaxies have a core radius 
re ~ C(300 pc), see Table IVnll We found that Tc is highly constrained if the stellar disk has 
a dominant behavior close to the center of the galaxy. However, we notice that even when re is 
different for each disk mass model there is an interval of values where Tc is consistent for all mass 
models within the la for most galaxies and some within the 2a error, as shown in Figures [2l [3l 
m andlH This provides us with enough tools to later test the Universe background evolution and 
large scale structure formation in order to continue testing the particle-physics-motivated BDM 



model. This energy transition is similar to the mean energy of a relativistic neutrino (see e.g. 3ll |) 
which is such that T,m^ < 0.58 eV(95% CL). 

The dispersion on Ec is much smaller than that of re and this is consistent with our BDM model 
since Ec is a new fundamental scale for DM model while re depends on r^, po which depend on the 
formation and initial conditions for each galaxy. We test the BDM profile with a fixed energy Ec 
finding that the proposed profile still being competitive compare with NFW profile. BDM profile 
has a rich -particle physics motivated- structure that, through its transition from CDM to HDM 
(when p ~ Pc), allows for a more physical explanation of the rotation curves of the different galaxies 
treated here. 
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We organize this Appendix as follow, in section[X]are shown the fitted values corresponding to the 
diet-Salpeter, for the NFW and BDM profile, and the free T^,, for BDM profile, mass models. In 
section [B] we show the fitted values for the minimal disk, min.disk+gas, Kroupa and diet-Salpeter 
when considering Burkert, Table IXTVl and pseudo-isothermal. Table [XVl profiles. Finally, in the 
section O we show the likelihood contour plots and the rotation curves for the different galaxies 
and mass models. 



A. DIET-SALPETER AND FREE MASS MODEL 



In this section we present the results of the diet-Salpeter mass model for the BDM, Eq. ([5]), 
and NFW, Eq. (Ilip . profile. We also present the results for the free T^ mass model for the BDM 
profile, as explained in the Section IVl 
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TABLE X: It shows the fitted values when considering the diet-Salpeter mass model for the value of the 
stellar disk. BDM and NFW parameters are show in columns (2-7) and (8-10) respectively. Units and the 

set of galaxies are as show in Table HIl 
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TABLE XI: Here we show the values obtained from the inner analysis for the diet-Salpeter mass model, 
with the profile pi„. All the parameters are explained in Table [V] The dashed (-) means that the Xred could 
not be computed since the number of data taken to do the analysis were too small. 
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TABLE XII: We show the main parameters, Vc, pc and Ec, for the BDM model as defined in Sec. Illll for the 
diet-Salpeter mass model. The values for G.A. galaxies are those obtained with the BDM profile, Eq. ([Sjl. 
The values for G.B. galaxies are those obtained from the inner analysis using pin as explained in Sec lVIl 
The cases with central value r^ < 10~® the energy value Ec cannot be computed, in this case we show a 

single dash. Units are as in Table. IVIII 
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^■»^-0.56 


8.95t?;?« 


0.41 


1.63 


NGC2903 


1 


r, or, + 0.02 

^"^^-0.03 


S.OOteia 


9 09+0.03 

^■■'^-o.os 


lo.ists:?! 


0.22 


1.29 


G-B NGC2976 


0.5 


. 9f,+0.06 


01+6.34 

«.Ul-6.33 


9 nq+0.07 


8.6lt"7:58 


0.03 


0.85 


NGC3198 


0.98 


c- 17+0.02 
''•-'-'-0.02 


S.OOlsl 


c ns+0.11 


10.17t«;«? 


0.21 


0.47 


NGC3521 


1 


1 01+0.02 

i-yj^_o.oi 


9.34t^7:« 


1 Q1 +0.005 
^•yi-0.005 


9.59tr.r8 


0.01 


1.1 


NGC925 


1 


Q nr,+0.08 


7.52tr.^^ 


lo.Mig:^^ 


< 1 


<0.01 


0.49 


NGC2403 





10 r, + 0.04 

-^-■^■^-0.04 


6.92+1^1 


0.0+0.001 


9 75+^-" 

y-'°-7.95 


0.19 


0.82 


G.C NGC5055 





7(5+0.01 
•^■'"-0.01 


,0+6.11 
».43_6.11 


0.0+0.01 


< 1 


<0.01 


1.26 


NGC7331 
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r, 19+0.001 

^•^^-0.001 


13 41+"-^ 


35.643l0.34 


10 24+0.1 

lU.Z1_g g-^ 


0.02 


0.43 



TABLE XIII: We show the fitted values when we let the mass-to-light ratio of the stellar disk as a free 
parameter. Columns (2-5) are the BDM parameters with units given in kpc for distances and Mo/kpc'' for 
the po- Column (6) is the logarithm of the fitted mass of the stellar disk and in (7) its associated mass-to- 
light. We can see that some galaxies are consistent with the minimal disk+gas, and must of them get values 
closer to stellar mass suggested by Kroupa instead of the Salpeter IMF. The set of galaxies are as shown in 

TableHI] 
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B. ISO AND BURKERT RESULTS 



In this section we present the tables where we show the results of the fit using the pseudo- 
isothermal (Eq. ()16p ) and Burkert profile (Eq. (jl4p ) for all the for mass models as explained in 
Section |Vl The values for BDM and NFW for min.disk, min.disk+gas and Kroupa mass model 
were presented in the Sec. |Vl] in Tables HU llllt llVi 



BURKERT 



Galaxy 
(1) 



MIN. DISK 



(2) 



logloPO Xred 

(3) (4) 



MIN. DISK + GAS 



'^rcd 



(5) 



logio PO 



(6) (7) 



KROUPA 
logio Po 



•^rcd 



(8) 



(9) (10) 



DIET-SALPETER 



(11) 



^rcd 



logio PO 
(12) (13) 



G.A 



DDO 154 
NGC 2841 
NGC 3031 
NGC 3621 
NGC 4736 
NGC 6946 
NGC 7793 



2.44 
3.23 



+0.02 



1.35 



+0.01 
-0.01 
+0.003 
0.003 
q qo+0.01 
■J-JJ-COl 
Q 29+0001 
'^■^''-0.001 
+0.01 



7.52 
8.95 



+5.62 



9.45 



+6.66 
-6.46 

+7.07 



7.08 

o 9R+5.96 
8-2b_5.gg 

+8.47 
47 



10.6^ 



2.86 
1.90 



+0.01 
-0.01 



8.63 

8.47 



+6.27 
-6.27 

+6.43 
-6.44 



0.57 

1.53 
4.34 
7.94 
1.73 
2.11 
6.11 



r, qi+0.02 

^•'-'^-0.02 
9 r,, +0.002 
^■^^-0.006 

1 r,r,+0.004 

^■^^-0.001 
n SQ+0003 

'-'■"^-o.ooi 

Q Q-+0.0001 

"■^^-0.0001 
.+0.002 



7 48+5-64 

1+7.42 

.55 



9.39"^ 



9.71 



+7.28 
-7.39 



1.77 



-0.006 



q 40+7.32 
^•^•-'-7.06 
j+10.84 
10.52 

+6.51 



1 77+0.01 
^•' '-0.01 



12.8] 
9.04 

8.49' 



+6.49 
-6.49 



0.53 
5.79 
5.45 
55.7 
15.8 
7.13 
5.68 



2 47+0.02 
^•^'-0.02 
+0.01 
0.01 



7A^+_ll^ 



4.08 
3.5__ 
5.77 
0.22 
8.12 
1 



+0.02 
0.03 
+0.02 
-0.02 
+0.05 

0.05 
+0.06 

0.06 
+0.01 
-0.02 



8.62 
8.22 
7.63 
10.5 
7.52 
8.35 



+6.34 
6.33 

+6.29 

6.3 

+5.51 

5.51 

+8.76 

76 

+5.52 

-5.52 

+6.44 



0.43 
3.02 
5.12 
5.59 
1.30 
1.14 
7.88 



9 c:o+0.03 
^■>J'-'_0.03 
Q ocr+0.03 
■^■°3_0.03 

1 90+°°^ 

Q 41+0.05 
^■^^-0.05 
+0.05 
0.05 



0.55"* 



309 



+322 
-105 



3.09 



+0.03 
-0.03 



7.41 
8.75 
9.06 
7.21 
9.80 
6.65 
7.98 



+5.59 

5.59 
+7.02 

7.01 
+7.82 

7.81 
+5.15 

5.17 
+7.40 

7.40 
+4.89 

-4.9 
+6.14 
-6.13 



0.39 

0.93 
4.37 
3.64 
3.33 
2.49 
4.30 



G.B 



IC 2574 
NGC 2366 
NGC 2903 
NGC 2976 
NGC 3198 
NGC 3521 
NGC 925 



,„ ^+0.19 
J^lJ-O_0.19 

1 02+0-04 
^-="^-0.03 
9 cc+O.Ol 
^■3O_0.01 
, 77+0.03 
^-' -0.03 



6.90 
7.81 
8.77 
8.40 



+5.06 

5.1 
+6.43 

6.3 
+6.59 

6.52 
+6.62 

6.63 



4.50 
2.18 

7.98 



+0.01 
0.01 

+0.01 
0.01 

+0.09 
-0.09 



1.0212 



+5.79 



5.83 
+6.87 



9-00I?;a 



7.35 



+5.55 
-5.47 



0.29 
1.63 
1.43 
0.95 
0.67 
1.25 
0.19 



8.71 
0.76 
2.43 
0.32 



+0.17 

0.17 

+0.01 

0.03 

0.01 

0.01 

+0.003 

0.003 



4 on+^-oi 
^•*J'J-o.oi 

2 14+0.01 
^•-^^-0.01 



6.86 
8.33 
8.81 
9.41 
8.09 



+ 5.13 

5.18 
+7.16 

6.82 
+6.60 

6.58 
+7.3 

07 
+5.90 

5.91 



9 01+'^-^** 

^■"^-6.98 



2.14 



+0.01 
-0.03 



.15 



+5.86 



0.31 
3.94 
1.31 
36.3 
0.42 
1.32 
13.7 



14.7 
1.75 
3.06_ 
>10'' 



+0.57 
0.54 

+0.04 
0.05 



+0.01 
0.01 



8.76 
24.5 
21.2 



+0.04 
0.04 

+ 1.32 
1.25 

+ 1.02 
0.95 



a K9+4.96 

7 79+6.37 
'•'^-6.37 

+6.37 

6.45 

+6.19 

-6.19 

7.271:^-22 

+5.14 
13 



8.55! 

7, 



6.77 



0.69 
1.34 
2.01 
1.18 
2.92 
8.64 
1.19 



20.8+i;iI 
1 70+0.05 

-■--'8-0.05 
9 ocr+O.Ol 

^■»a_o.oi 

+3.06 
3.06 
q 74+0.05 
y-'^-0.05 
+ 135 
50.1 
+ 18.5 
-8.61 



23. 7"^ 



162] 
61.1 



6.51 
7 

8.63 
7.50 



+4.89 

4.89 

+6.36 



-6.35 
,+6.44 
j.43 
1+5.39 
-5.39 
7 1 7+5.16 
'-J^'-5.13 
f, 99+5.02 
"■^^-5.01 

6-54ttI^3 



1.16 
1.34 
1.23 
5.29 
0.93 
11.5 
2.54 



NGC 2403 

G.C NGC 5055 

NGC 7331 



2.78 
2.26 
2.36 



+0.01 
0.01 

+0.01 
0.01 

+0.01 
-0.01 



8.32 
8.91 
9.05 



+5.83 

5.91 

+6.59 

6.58 

+6.88 

-6.88 



4.33 
6.34 

1.47 



0.84 
0.07 
2.16 



+0.001 

0.001 

+0.0001 

0.0001 

0.01 

0.01 



9.42_ 
12. 5J 
9.12 



+7.45 



+6.96 



27.8 

94.4 
1.08 



3.S 
23 
17.4 



+0.01 
-0.01 

+0.27 

-0.26 

0.19 

0.19 



7.98 



7.23 



+5.63 
5.63 
+4.79 
-4.8 
+5.33 
-5.34 



1.6 
4.11 

7.06 



3.9 
2.52 



+0.01 
0.01 

+0.03 
0.03 



27.58ir4l 



7.96 
8.79 
7.00 



+5.59 
5.62 

+6.89 
6.89 

+5.15 
5.14 



2.25 

3.47 

15.57 



TABLE XIV: We present the fitted values of tfie free Burkert parameters, Eq. (|14|l . for all the mass model. 
Each column. The r^ is given in kpc. Logarithm base 10 of the central densities po is given in MQ/kpc^. 
The value of x^ is normalized to the numbers of data points minus the number of free parameter. The set 

of galajcies are as show in Tablellll 



C. ROTATION CURVES AND LIKELIHOODS CONTOURS 



In this Sec. we present the figures where we plot all the different profiles, presented in the mass 
model Sec. IV D) the confidence level contours obtained from the numerical fitting as well as the 
consideration made for each galaxy. The groups are described in the Sec. I VII In each subsection 
we present the rotation curve for each galaxy and mass model, at the end of each subsection we 
show the respectively confidence level contours. The G.A. galaxies are presented in Sec. IC C.Al in 
subsection Sec. ICC.BI and IC Al are presented the galaxies included in G.B for the BDM and the 
BDM-inner analysis. Finally, in subsection IC C.Cl we present G.C galaxies. 
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PSEUDO-ISOTHERMAL 



Galaxy 
(1) 



(2) 



MIN.DISK 
logpo 
(3) 



MIN.DISK+GAS 



(4) 



y2 



(5) 



logpo 



(6) (7) 



KROUPA 

logpo 



DIET-SALPETER 



y2 

'^rod 



(8) 



(9) (10) 



(11) 



'^red 



logpo 



(12) (13) 



DDO 154 
NGC 2841 
NGC 3031 
G-A NGC3621 
NGC 4736 
NGC 6946 
NGC 7793 



1.30'' 



+0.01 
0.01 

0-^4+0.01 
~ 10 



7 c:c + 5.6 



0.01 
7 



1.01 



+0.01 
-0.01 



2 63+°°^ 

1 oa+o-''^ 



0.95 



+0.01 
0.01 



11.0 

22.80 

8.65 

8.10 

8.88 
;.52 



+8.61 

8.53 
+ 20.4 

20.5 
+6.37 

6.35 
+6.10 

5.99 

+6.52 

-6.52 

+6.48 

-6.49 



0.50 
5.36 
6.98 
0.86 
168 
1.66 
5.74 



1.21^ 



0.01 
0.01 

0.82lHi 



7 trr,+ 5.68 

'•"'^-5.69 
q 47+7.32 

,+7.14 
6.26 



0.58 
0.38 
1.68 
0.86 



+0.01 
-0.01 
+0.01 

0.01 
+0.06 

0.03 
+0.01 

0.01 



9.08 
9.05 
9.51 
8.53 
8.56 



7.15 
6.43 
7.64 
7.22 

+6.27 
6.06 

+6.56 
-6.56 



0.40 
10.9 
25.6 
1.89 
47.2 
4.09 
5.22 



^■■'^-o.oi 

< 0.01 

+0.01 



1.42 
1.99 



+0.01 



< 0.01 



4.79! 
0.87 



+0.01 
-0.01 



5.64 
5.63 
11.3 
11.3 

+6.50 

6.49 

+5.83 

-5.82 

21 0+^^-^ 

7 40+5.49 

+6.54 
6.55 



7.45- 

13.6 

8.42' 
7.94 



8.44 



0.36 
0.72 
5.14 
2.8 
5.96 
1.10 
7.02 



1 36+001 

1-3D_0.01 

< 0.01 

>108 



7 4Q+5.61 
'•^'^-5.61 

14.95+12-4 



+0.03 



4.56 
< 0.01 



41.6 
1 



+ 7.14 

4.87 

+0.02 

-0.02 



6.54 
7.29 
12.0 
6.66 
7.94 



+4.79 

-5.76 
+5.23 



+ 10.7 



+4.91 
-4.91 
+6.09 



0.35 
3.14 
20.4 
3.01 
4.29 
2.48 
4.78 



G.B 



IC 2574 
NGC 2366 
NGC 2903 
NGC 2976 
NGC 3198 
NGC 3521 
NGC 925 



6.62 
1.16 



+0.12 
0.11 

+0.02 
0.02 



og+o-oi 
^■^^-0.02 

-I -19+0.02 
-^■-^^-0.02 
+0.02 
0.02 
+0.01 
0.01 
+0.06 
0.06 



1.37; 

0.90 
5.01 



6.84 
7.75 
9.96 
8.35 
8.40 
9.16 



+5.01 

5.04 
+6.24 

6.23 
+ 7.73 

7.77 
+6.56 

6.57 
+6.17 

6.22 

7.12 

7.12 



7 29+5.49 



0.26 
2.07 
6.47 
1.03 
2.79 
5.70 
0.20 



5.47 
0.06 
0.51 
0.04 
0.36 
0.63 
1.03 



+0.11 
0.11 

+0.01 
0.02 

+0.02 
0.02 

+0.01 
0.01 
0.01 
0.01 

+0.04 
0.02 
0.01 
0.01 



6.81 
9.35 
9.46 
10.4 
9.47 
9.43 



+ 5.07 
5.12 

+8.04 
8.35 
7.35 
7.13 

+9.08 
.36 

+6.97 
7.58 

+7.03 
7.76 



8.20l6;g 



0.27 
11.8 
7.58 
36.7 
8.04 
7.19 
10.3 



8.95 
1.36 
0.59 



+0.30 
0.31 

+0.02 
0.02 

+0.02 
-0.02 



> 10" 



6.57 
7.47 
9.31 
7.69 



+4.91 
4.9 

+5.83 
6.45 

+ 7.13 

-7.12 
5.61 
5.56 



3.80 

15.38 

12.4 



+0.02 
0.02 

+0.83 
0.79 
+0.54 
0.50 



+5.34 
5.34 



6.6"' 



6 73+5-10 
"•'■^-5.09 



0.64 
1.79 
6.04 
1.15 
3.42 
8.62 
1.13 



11.7 
1.06 
0.56 
3.00 
6.70: 
46.3 
> 10 



0.64 
0.57 

+0.03 
0.03 

+0.02 
0.02 

+0.01 
0.01 

+0.04 
0.04 
11.7 
7.39 

6 



6.48 
7.62 

9.35 
7.56 
6.96 
6.24 
6.47 



+4.86 



+6.34 
-6.29 
+ 7.16 
7.24 
+6.20 
15 



+4.68 
-4.98 
+5.03 
-5.03 
+0.60 
-0.60 



1.12 

1.64 
6.46 
6.01 
6.40 
11.49 
2.47 



NGC 2403 

G.C NGC 5055 

NGC 7331 



o.84t°:gl 

0.14«-.oi 

+0.01 
0.01 



0.35 



8.72 
10.6 
10.0 



+6.23 
6.31 

+8.24 
8.24 

+ 7.91 
7.73 



1.04 
2.79 
1.80 



0.57 



+0.01 
-0.01 



< 0.002 14.2 



9 00+'^-*" 
y-uu_g 3g 

+11.9 
11.9 



0.20 



+0.01 



10.5 



+8.35 

-8.23 



1.14 

3.43 
2.18 



1.53 



+0.01 
-0.01 



15 04+0.18 
J^3-U^_0.18 

> 10^ 



8.19 
6.53 
6.73 



+5.85 
5.84 

+4.72 
4.72 

+4.66 
5.02 



0.85 
4.14 
10.2 



1.91 
45.3 

>108 



+0.01 
-0.01 
-1.69 

-1.58 



8.01 
5.99 
6.64 



+5.67 
5.68 

+4.29 
4.29 

+4.98 
4.70 



0.8 
13.3 
25.9 



TABLE XV: We present the fitted values of the free pseudo-isothermal parameters, Eq. (|16p . for all the 
mass model. Each column. The re is given in kpc. Logarithm base 10 of the central densities po is given in 
Mq /kpc^. The value of x^ is normalized to the numbers of data points minus the number of free parameter. 

The set of galaxies are as show in Table HTl 



A. Inner analysis 

We present the different G.B. galaxies that where analyzed when the DM halo is the only mass 
component with the inner approximation for the BDM profile pin = 2pc(l + r/rc)~^ up to a radius 
r < Tjj with a constant positive slope shown in Fig. ([5]-[8|) with a purple, long, dashed line. From 
this first approach we obtain the the core distance Vc and the central density pc, then we use these 
parameters to obtain po and r^ from the fit of BDM profile to all the set of observational data 
(thick black line). Fitted values are presented in Table. |V] and the confidence levels between pc 
and Tc in Fig. [3l 



C.A. GROUP A 



In this section we present the considerations made in each one of the galaxies for the numerical 
analysis, this galaxy has in particular fitted values of Tc < r^ when analyzed with the minimal disk 
model, we have named this set of galaxies as G.A. The conclusion are presented in Sec. IVIII At 
the end of the section can be seen the confidence level for each galaxy in every mass model. 
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Hin. Disk 



Min. Disk+gas 



Kroupa 



diet-Salpeter 




6.92 
6.86 




0.62 1.19 1.77 2.34 



8.6016 

0.00 0.04 0.08 0.12 



0.19 0.54 0.89 1.24 




8.6687 
8.6526 
8.6365 

0.00 0.04 0.07 0.11 



10.11 
10.08 
10.05 
10.03 



0.05 0.15 0.24 0.34 



0.64 0.83 1.01 1.20 





7.4929 

0.00 0.03 0.05 0.0 



0.00 0.05 0.09 0.13 




0.01 0.04 0.07 0, 



0.34 0.43 0.53 0.63 



8.0048 



8.0916 

8.0744 



0.02 0.10 0.17 0.25 



0.04 0.12 0.20 0.28 




0.01 0.05 0.09 0.13 0.01 0.05 0.09 0.13 



7.05 
6.97 




0.32 0.87 1.41 1.95 




S.1614 

0.00 0.03 0.06 0.09 




0.00 0.11 0.22 0.33 




6.6920 

0.0001 0.0162 0.032. 




0.02 0.13 0.24 0.35 




0.28 0.96 1.63 2.31 




7.1272 

0.0001 0.0122 0.0244 




0.48 1.22 1.96 2.71 



8.0520 
8.0374 
8.0227 

0.00 0.03 0.05 0.07 



<10' 7.25x10° 



0.00 0.02 0.04 0.06 




3.5314 3.5317 3.5319 



Tc [ kpc ] 

FIG. 2: We present the 2D likelihood contours for the BDM parameters po and re for galaxies with Tc ^Ts, 
group A. The two regions in each figure correspond to la (68%) and 2a (95%) confidence level. Columns 
from left to right correspond to minimal disk, minimal disk + gas, Kroupa and diet-Salpeter mass models. 
This clearly shows how the stars has an important contribution close to the center of the galaxy and erase 
trace of the core, we can see that even if the value of Vc that minimize x^ is zero, is consistent with values 

different from zero up to la. 



Min. Disk 



Min.Disk+gas 



Kroupa 




6.22 19.92 33.61 47.31 




!.12 3.15 4.13 




1.45 1.77 




2.38 3.25 4.12 




2.96 3.53 4.10 4.67 




2.00222 2.00223 2.00224 




6.29 11.59 16.8 




4.41 14.40 24.40 34.39 




2.89 3.90 




1.45 1.70 1.94 




0.55 1.49 2.42 




3.31 3.81 




L29 2.74 3.19 




7.30 19.64 31.S 




S.64 25.98 48.3 




2.34498 2.34498 2.34498 




79 1058.00 2109.00 3159.00 





1.39 13.75 26.11 38.47 




100.90 194.70 288.50 
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FIG. 3: This table shows the 2D hkelihood contours for the BDM parameters po and re for galaxies with 
re = Ts, group B. The different colors in each figure represent the la and 2<j confidence levels. Columns from 
left to right corresponds to minimal disk, minimal disk + gas, and Kroupa. When more mass components 

nrti inr*liirltirl in fVit> nnnKreiQ fVio f>r»n^r^iire r»f fVit> nnnfirl*:inr'ei lenT-tilG Kcmnmo Inrcrtir rticnlfincr in ^\\e^ Hiffinnl+ir 



33 



o 

M 

Q 



Ci 



o 



7.520 



7.514 



7.508 



0.000 



8.389 



Min. Disk 



0.012 



0.024 



0.00 0.04 0.09 0.13 



Min. Disk+gas 





7.575 



7.569 



7.56; 




0.000 



8.441 



8.434 



8.428 - 



0.000 



0.013 



0.007 



0.025 




0.014 




Kroupa 



4.144 




0.002 0.004 




16.53 24.78 33.03 



2.19 



2.14 



2.09 



2.04- 




1.15 3.21 5.27 7.33 



re [kpc] 

FIG. 4: This table shows the 2D hkehhood contours between the BDM parameters po and re, corresponding 
to 1(7 and 2(7 for galaxies with re = 0, group C. Columns from left to right correspond to minimal disk, 
min. disk+gas, Kroupa, and diet-Salpeter mass models. 
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FIG. 5: IC 2574 and NGC 2366 inner analysis. The gray points are the observations with its error 
bars. Below we show the difference between the observational and the theoretical approach, the 
purple region represent the error bars and the line the fitted curve. 
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FIG. 6: NGC2976 inner analysis 
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FIG. 7: NGC3198 inner analysis 
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FIG. 9: The 2D likelihood contours for the BDM parameters po and r^ for galaxies with re = rs, group B. 
The different colors in each figure represent the la and 2g confidence levels. Columns fi-om left to right 
corresponds to minimal disk, minimal disk + gas, and Kroupa mass models obtained with the inner analysis. 
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FIG. 10: Graphical visualization of the rotation curves for the galaxies DDO 154 belonging to the group G. A. 
where r^ < rs for all mass models. Is a dwarf galaxy, one of the first galaxies used to illustrate the conflict 
with the theoretical prediction of the ACDM [17] . No significant color gradient in the colors (B- V) and (B-R) 
is detected, constant color as function of radius in our models is well assumed with the value for T* = 0.32. 
For a single exponential disk we use value of fio = 20.8 mag arcsec? and Rd — 0.72 kpc. BDM is the best 
option since it predict a value of T*. very close to the minimal disk and gives the best fit. The values for 
Ts and re are very reasonable except when we considered the Kroupa IMF. From left to right and top to 
bottom, images from the fit considering minimal disk, minimal disk+gas, Kroupa, diet-Salpeter. In the top 
images of each galaxy are represent the four DM profiles: BDM, red squares; NFW, blue circles; Burkert, 
green diamond; Isothermal, black triangles. The dotted line is the gas component, the short dashed line is 
the stellar contribution and the tick points are the observational data with its respectively error bars. At 
the bottom images we plot the difference between the observed and predicted curves where the purple region 
represents the error of the observations and each different line represent the fitted curve of the different 
profiles with the observations taking into account already all the mass contributions. 



37 



R [ITJC] 




(a) Minimal disk 



(b) Min. disk + Gas 



R [I-pc] 



R [ITJC] 




30 40 



R[kpc] 



(c) Kroupa 



R[kpc] 



(d) diet-Salpeter 



FIG. 11: We present the rotation curves for the galaxy NGC 2841. Colors and symbols are as in Fig llOl This 
galaxy shows a inner component of characteristic radius of Rd — 0.72kpc but we consider it as one-component 
with total mass of the star disk of A'h = lOn.O4M0 and Ra = 0.72kpc. BDM and NFW presents the best 
fitting when considering the DM as the only one which contributes to the rotation curve. BDM does not 
differ much from NFW since Vc is very close to zero. Fits when we considered other mass elements tends to 
zero the value of the core, but are consistent with la confidence intervals. From left to right, images from 
the fit of the galaxy NGC 2841 considering minimal disk, minimal disk+gas, Kroupa, diet-Salpeter. 
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FIG. 12: We display the rotation curves for the galaxy NGC 3031. Colors and symbols are as in Fig llOl This 
galaxy is the proto- typical grand design spiral. The profiles show evidence for a central small component. 
To describe the outer disk we use the parameters /^o = 15.9 and Rd — 1.94. The (J - K) color profile shows 
no significant color gradient. We assume a constant T^, = 0.8 for the disk. We can not extract information 
of these galaxy since we do not have observational points close to the center. From left to right, images from 
the fit of the galaxy NGC 3031 considering minimal disk, minimal disk+gas, Kroupa, diet-Salpeter. 
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FIG. 13: The rotation curves for the galaxy NGC 3621. Colors and symbols are as in Fig llOl There is 
no a convincible evidence to believe that this galaxy has a bulge, neither the color gradient, morphology 
or deviation for a single exponential disk. Therefore we assume a single-component disk to larger radii 
using an exponential of mass M* = LOlO^^Afo and Rd — 2.61 kpc. It is worth noticing how close are 
the predicted values of the BDM and NFW to the expected values for the maximum stellar disk and at 
the same time they provided the best fit to this galaxy. The value of the core is underestimated the more 
mass components are considered. Since the galaxy present sufficient data we proceed with the inner analysis 
obtaining consistent results with the minimal disk analysis. From left to right, images from the fit of the 
galaxy NGC 3621 considering minimal disk, minimal disk+gas, Kroupa, diet-Salpeter. 
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FIG. 14: The rotation curves for the galaxy NGC 4736. Colors and symbols are as in Fig llOl The stellar 
component of the galaxy has a complex structure, we analyzed it as one and two component stellar galaxy, 
where its central component has a the parameters /lo ~ 11.8magarcsec~^ y h = 0.26kpc. For the outer disk 
we assume a mass of M* = 1.8 * lO'^" Mq and Rd = 1.99 kpc. The galaxy exhibits a small color gradient 
and we assume a constant T* = 0.6 and T* = 0.9 for the outer and inner disk respectively. The best fit is 
given by the minimal analysis and is consistent with Kroupa IMF without bulge. The analysis with bulge 
does not improve significantly the value of x^ ■ From left to right, images from the fit of the galaxy NGC 
4736 considering minimal disk, minimal disk+gas, Kroupa, diet-Salpeter. 
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FIG. 15: Here we display the rotation curves for the galaxy NGC 6946. Colors and symbols are as in 
Fig llOl The galaxy present a small central object but does not improve remarkably the analysis, so we treat 
the stellar disk as one component using a mass of M* = 5.8 * IO^OAIq and Rd ~ 2.97kpc. From left to 
right, images from the fit of the galaxy NGC 6946 considering minimal disk, minimal disk-|-gas, Kroupa, 

diet-Salpeter. 



42 




(a) Minimal disk 



(b) Min. disk + Gas 



R [I-pc] 



L 
120 

m 


12 3 4 5 6 7 

^,,,1111"""! IHl!lii)^ii^ 


so 


jjiJ^^°^^ 


60 


ii^^ 




i'L^ ^ -- ^ ' \ ^ X 


40 


W -■' ^' --•- ^''~.,' - 


20 


if ' , ■■■■■■■■■■■ 




R[I-pc] 
12 3 4 5 6 7 



i„iiii"""ililiiiiiii[ii!Jni]j 




(c) Kroupa 



(d) diet-Salpeter 



FIG. 16: Here we display the rotation curves for the galaxy NGC 7793. Colors and symbols are as in Fig llOl 
This galaxy is a late-type Sd spiral galaxy. It shows a nuclear star cluster, nevertheless its photometric and 
dynamical importance has become negligible. The (J-K) profile has a small color gradient, so we proceed 
taking an average of the T* for the inner and outer parts of the stellar disk and assume a constant value for 
T* = 0.31. All the profiles roughly predict equal results, close to Kroupa's value, even so, BDM offers the 
best fit for these galaxy with physical plausible parameters values. From left to right, images from the fit of 
the galaxy NGC 7793 considering minimal disk, minimal disk+gas, Kroupa, diet-Salpeter. 
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C.B. GROUP B 



In this section we present the considerations made the galaxies with fitted values of Tc = r^ when 
analyzed with the minimal disk model, we have named this set of galaxies as G.B. The conclusion 
are presented in Sec. IVIIi At the end of the section can be see the confidence level for each galaxy 
in every mass model. 
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FIG. 17: The rotation curves for the galaxy IC 2574. Colors and symbols are as in Fig llOl A dwarf galaxy. 
We analyze this galaxy as a one component exponential disk with mass M* = 1.04 * 10 A/© and Rd = 2.56 
kpc. We assume a constant mass-to-light ratio with value T* — 0.44. All the profiles predict the minimal 
plausible value for T^. All the cored profile adjust equally well and NFW fit specially bad. In these kind 
of galaxy we have the best fit for r^ = r^ even so we do an analysis using the inner behavior of the BDM 
profile. These galaxy has a very smooth slope and we can not distinguish where the maximum could it be, 
so we take all the points and adjust it to the inner BDM profile in order to find a plausible value for the 
core and the central density re — 3.2kpc and pc = 3.4 * lO^MQ/kpc'^ for the minimal analysis. Since we take 
all the points to do de inner analysis it is very difficult to break the degeneracy between the parameters po 

and Tc. 
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FIG. 18: We display the rotation curves for the galaxy NGC 2366. Colors and symbols are as in Fig llUI 
For these galaxy is clear that non circular velocities are affecting the rotation curve. We do inner analysis 
taken data below the 2.2 kpc obtaining values for the core Tc = 7.4kpc and the central density pc = 
2.4 * lO^MQ/kpc^. even if the parameters are not absolutely correct can be taken into account in order 
to do predictions of the behavior of the curve. For the stellar disk we use a, Rd = 1.76 kpc and mass 

M* = 2.5*1O*M0. 
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FIG. 19: We present the rotation curves for the galaxy NGC 2903. Colors and symbols are as in Fig llOl The 

profile of this galaxy can be decomposed into two components, but the central object provides a small effect 

on the rotation curve fit and it is ignored. For the outer disk we take values of M^ = 3.1 * W^OMq and 

Rd = 2.40 kpc. The inner analysis is carried out taking data points below 4.5 kpc and obtaining values for 

core and central density of re = 1.7kpc and pc — 6.7 * lO^Afo/kpc'^. 
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FIG. 20: We display the rotation curves for the galaxy NGC 2976. Colors and symbols are as in Fig llOl The 
galaxy show a color gradient, for sake of simplicity we take T, = 0.52 for all radii. We assume the value 
Ho = 17.78 mag arcsec"^ and Rd — 0.91 kpc. For the inner analysis we take all the observational data, 
up to 2.5 kpc, since the slope of the rotation curve is very smooth and can not be distinguished where the 
maximum velocity could be located. The fitted values obtained from the inner analysis are re — 0.25kpc 
and pc = 1.3 * lO^Mo/kpc'^. It is difficult to break the degeneracy between re and po since we have no more 

data. 
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FIG. 21: We present the rotation curves for the galaxy NGC 3198. Colors and symbols are as in Fig llOl This 

galaxy shows a complicated stellar disk in the center, it could be assume two component but for simplicity 

and because it provides a small effect on the rotation curve we analyze it with only one component with 

parameters, /lo = 17. lmagarcsec~^ and Rd = 3.18kpc. For the inner analysis we take data below 6.8 kpc, 

and obtained values of r^ = 0.21kpc, r^ = 16.64 and po = 7.7 * lO^MQ/kpc^ for the minimal disk. 
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FIG. 22: The rotation curves for the galaxy NGC 3521. Colors and symbols are as in Fig llOl Because the 
nature of the galaxy it is difficult to obtain reliable data from the stellar component (see deBlokOS). even 
so, we treat the stellar disk as a single component with mass A/* = 1.2 * IO^^Mq and Rd = 3kpc. A color 
gradient is present. The inner analysis is carried out with data below the 3.7 kpc and obtaining values for 
the core and central density Vc = 2.4 kpc and r^ — 5.5 kpc and po ~ 2.4* lO^Mo/kpc^ for the minimal disk. 



49 




« -5 



^ \^i:"::= s :::i»''' 



R[kpc] 
(a) Minimal disk 

R[I-pc] 



R[kpc] 

(b) Min. disk + Gas 

R[I-pc] 




R[kpc; 
(c) Kroupa 



R[kpc] 
(d) dict-Salpeter 



FIG. 23: We display the rotation curves for the galaxy NGC 925. Colors and symbols are as in Fig llOl It 
is classified as a late- type barred spiral. Shows no evidence for a bright central component. There is also 
no evidence for a strong color gradient and we assume a T* = 0.65 for the stellar disk. The stellar disk is 
described by a central surface brightness fio ~ 18.9 magarcsec" and a scale length Rd = 3.31 kpc. This 
galaxy is specially difficult to fit with reasonable values specially because it present a very smooth slope. 
For inner analysis we take data below 6 kpc, obtaining values for the minimal disk analysis of Vc = lO.Skpc 
for the core and pc = 2 * lO^^A^o/kpc^ for the central density. Unlike the previous IC 2574 and NGC 2976 
galaxies we have few data points after the 6 kpc so we can break the degeneracy and obtain values for 
Ts = 24.6 kpc and po = 8.1 * lO^MQ/kpc^ for the minimal disk. 
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C.C. GROUP C 

In this section we present the considerations made for the galaxies with fitted values of Tc = 
when analyzed with the min.disk mass model. We have named this set of galaxies as G.C . The 
conclusion are presented in Sec. IVIII for further discussion. At the end of the section can be see 
the confidence level for each galaxy in every mass model. 
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FIG. 24: The rotation curves for the galaxy NGC 2403. Colors and symbols are as in Fig llOl It is a late-type 
Sc spiral. The galaxy shows a (J-k) color gradient that produce a small difference in terms of T* for the 
inner and outer part of the stellar component. Since the color gradient is small we considerer a constant 
T* = 0.41. The surface brightness profile shows evidence for a second component in the inner parts but 
it is uncertain if this constitutes a dynamically separate component. For our purpose we consider a single- 
component in the stellar disk described by a central surface brightness fio = 16.90mag arcsec"^ and a scale 
length Rd = 1.81 kpc. For the number of inner data points less than 2.1 kpc it is clearly that we can do 
a inner analysis as for B group galaxies, in which the results are re = 0.004kpc and pc = 2.1 * 10* wich is 
consistant until la for the confidence level with zero. From left to right, images from the fit of the galaxy 
NGC 2403 considering minimal disk, minimal disk+gas, Kroupa, diet-Salpeter. 



51 



RPqjc] 



R [Iqjc] 




R[kpc] 

(a) Minimal disk 

R [Iqjc] 



R[kpc] 

(b) Min. disk + Gas 

R [Iqjc] 







10 


20 


30 


40 






*H A^ 








^•« 


-50 


w> 


^•■%»^ 


h/*"''' 


■^ "•*'■■•■ 




-100 


V 










-ISO 













10 20 30 

R[lq)c] 




(c) Kroupa 



(d) diet-Salpeter 



FIG. 25: Here we display the rotation curves for the galaxy NGC 5055. Colors and symbols are as in Fig llOl 
We take the parameter for this galaxy as /xo = 16.7 mag arcsec"^ and Rd = 3.68 kpc for the stellar disk. 
Inner analysis is in agreement with the zero value for the core until la of confidence level. From left to 
right, images from the fit of the galaxy NGC 5055 considering minimal disk, minimal disk+gas, Kroupa, 

diet-Salpeter. Colors and symbols are as in Fig llOl 
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FIG. 26: The rotation curves for the galaxy NGC 7331. Colors and symbols are as in Fig llOl We assume 
values for the parameters /xo = 15.3 mag arcsec"^ and Rd = 2.42 kpc for the outer disk, respectively. We can 
not extract relevant information from these galaxy because the observational data is not sufficiently close to 
center. From left to right, fits considering minimal disk, minimal disk+gas, Kroupa, diet-Salpeter. Colors 

and symbols are as in Fig llOl 
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